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Abstract 

Attempts to construct chromodyons — objects with both magnetic charge and non-Abelian 
electric charge — in the context of spontaneously broken gauge theories have been thwarted in the 
past by topological obstructions to globally defining the unbroken non-Abelian "color" subgroup. 
In this paper we consider the possibility of chromodyons in a theory with SO (5) broken to SU(2) x 
U(l), where the topological obstructions are absent. We start by constructing a monopole with 
only magnetic charge. By exciting a global gauge zero mode about this monopole, we obtain 
a chromodyonic configuration that is an approximate solution of the field equations. We then 
numerically simulate the time evolution of this initial state, to see if it settles down in a stationary 
solution. Instead, we find that chromo-electric charge is continually radiated away, with every 
indication that this process will continue until this charge has been completely lost. We argue that 
this presents strong evidence against the existence of stable chromodyons. 
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I. INTRODUCTION 



The magnetic monopoles that arise in spontaneously broken gauge theories can easily be 
generalized to dyons that have a U(l) electric charge in addition to their magnetic charge. 
It then is natural to ask whether, in cases where the unbroken symmetry is non-Abelian, it 
is possible to have monopoles carrying non-Abelian electric charge. Such objects, referred 
to as chromodyons, were first considered in the context of an SU(5) grand unified theory. 
After attempts to construct such chromodyons failed [1] , it was shown that the non-Abelian 
magnetic charge of the SU(5) monopole creates a topological obstruction to the existence of 
non-Abelian "color" electric charge [2-7], and the issue was abandoned for a number of years. 
Since then, however, it has been realized that, with other choices of gauge group, there can 
be monopoles with purely Abelian magnetic charge, even though the unbroken gauge group 
is non-Abelian [8]. Because there is then no topological barrier to a color electric charge, it 
is natural to revisit the subject and consider whether chromodyons can exist. 

Let us first recall how ordinary U(l) dyons arise. When there is an unbroken U(l) 
symmetry, any soliton with a nonvanishing charged field has a zero mode corresponding to 
a shift in the phases of the complex charged fields. Exciting this mode in a time- dependent 
fashion produces a U(l) charge. If the U(l) symmetry is gauged, it may be possible to 
gauge transform away the time-dependence of the phase, but the gauge-invariant electric 
charge remains. The simplest example of this occurs with SU(2) broken to U(l), where the 
Julia- Zee dyon [9] arises from rotation of the phase of the massive vector boson fields in the 
core of the 't Hooft-Polyakov monopole [10, 11]. 

Similarly, if there is an unbroken non-Abelian symmetry, excitation of the gauge orienta- 
tion zero modes of a soliton gives rise to a non-Abelian electric charge. In the theory with 
SU(5) broken to SU(3) color x U(1) EM , the unit monopoles have nontrivial fields that are not 
invariant under the unbroken SU(3). Hence, one would expect to be able to generate dyons 
that were charged under the color SU(3) (hence the term "chromodyon" ) by exciting the 
resulting global gauge zero modes. However, Abouelsaood [1] found that, because the gauge 
potential has a 1/r tail in the unbroken subgroup, some of the expected zero modes are non- 
normalizable, and the proposed construction does not go through. A deeper explanation for 
this was given by Nelson and Manohar [2], and by Balachandran et al. [3-5], who showed 
that the non-Abelian Coulomb magnetic field creates a topological obstruction that prevents 
one from globally defining a basis for the unbroken color subgroup. This inability to define 
"global color" is the fundamental reason for the nonexistence of the SU(5) chromodyons 1 . 

These barriers to the existence of a chromodyon would both be absent if the total magnetic 
charge were purely Abelian. This can certainly be achieved by assembling a collection of 
magnetic monopoles such that the non-Abelian components of their charges sum to zero; a 
two-monopole example of this was studied by Coleman and Nelson [13]. However, what we 
need to produce a chromodyon is a single monopole with purely Abelian magnetic charge. 
While there are no such monopoles in the SU(5) theory, they do exist in a theory with 
SO(5) broken to U(l)xSU(2). These were first discovered [8] in the BPS limit, where they 



The SU(5) monopoles do have dyonic counterparts with an electric charge in the U(l) subgroup defined by 
the magnetic charge, which lies partly in the unbroken SU(3) [12]. However, because the electric charge is 
restricted to this subgroup, implying that one cannot generate full color multiplcts of states, and because 
the color electric charge is strictly proportional to the Abelian electric charge, these are chromodyons only 
in a limited sense. 
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appear as spherically symmetric classical solutions that are characterized by a "cloud size" 
b that can take on any positive value. These can be interpreted as being composed of a 
massive monopole, carrying both Abelian and non-Abelian magnetic charge, and a massless 
monopole, with only non-Abelian magnetic charge. At the semiclassical level, the latter is 
manifested as a cloud of non-Abelian field, of radius b, that surrounds the core of the massive 
monopole and completely shields the non-Abelian part of its magnetic charge. In the BPS 
limit the energy is independent of b. However, if the Lagrangian includes a nonvanishing 
potential, the cloud size is no longer arbitrary, but rather is fixed. This then gives a magnetic 
monopole whose long-range field lies only in the U(l) sector, but whose core transforms under 
the unbroken SU(2) and thus gives rise to the gauge zero modes from which we might hope 
to construct a chromodyon. It is this system that we will study. 

We start, in Sec. II, by constructing the static SO(5) monopole. An analytic solution 
exists for the BPS case. However, as we explain later on, the possibility of varying the 
cloud size makes these unsuitable for our purposes. Instead, we must take for our starting 
point a non-BPS monopole, for which the field equations must be solved numerically. Then, 
in Sec. Ill, we construct a chromodyonic configuration from this monopole by applying an 
SU(2) gauge rotation and solving for the A field that is required by Gauss's law. Although 
this is not an exact static solution of the field equations, one might expect it to be close to 
the desired chromodyon. We test this by numerically simulating the time evolution with this 
as the initial configuration. We describe the details of this simulation in Sec. IV. The results 
are described in Sec. V. We find that, rather than evolving toward a stable chromodyon, the 
chromodyonic configuration continually radiates non-Abelian charge. Although we are not 
able to continue the simulation long enough to verify that this charge is completely radiated 
away, every evidence indicates that this will be the case. Our conclusions are summarized 
in Sec. VI. There are two Appendices containing some technical details. 

II. SO(5) MONOPOLES 

We are interested in theories with Lagrangian densities of the form 

C = — -Tr F^F^ Tr _D M <E> — V($) . (2.1) 

Here the gauge field and the adjoint representation Higgs field <£> are both written as 
imaginary antisymmetric 5x5 matrices. 

To describe the components of these and other adjoint representation fields, we will adopt 
the following conventions. In the defining representation, the generators of SO (5) are the 
ten 5x5 matrices 

( J mn )ij = -i(S im 5 jn - 5 in 5 jrn ) , 1 < m < n < 5 . (2.2) 
From these we can define six matrices 

K = \ {^e abc J bc + J a4 ) , k a = i Qe a6c J bc - J tt4 ) (2.3) 

that generate SO(4)=SU(2) xSU(2). We can then decompose any adjoint representation 
field P in terms of two triplets, Pfa and Pfy, and P^ (/i = 1, 2, 3, 4) via 

P = P (1) .h + P (2) -k + Pg ) ^ 5 . (2.4) 
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We will refer to P(i), P( 2 ), and P( 3 ) as the first-, second-, and third-sector components, 
respectively. 

SO(5) can be broken to SU(2)xU(l) in two inequivalent ways. In the first, corresponding 
to the decomposition SO (5) D SO (3) x SO (2), the SU(2) is the subgroup, with the generators 
J ah — e abc{h c + k c ) , that rotates the first three components of a five- vector among themselves. 
We will be concerned with the second possibility, in which the unbroken SU(2) is one of the 
factors of the SO (4) = SU(2) x SU(2) subgroup that mixes the first four components of a 
five- vector among themselves; we will choose it to be the subgroup generated by the k a . 

We will be seeking spherically symmetric monopole solutions. If we also require that the 
fields have positive parity, the most general spherically symmetric ansatz can be written as 2 

= £ai m r m A{r) , = f a H{r) , 

A a i{2) = e aim f m G(r) , $^ 2) = f a K{r) , 

A?, 3) = V2 F(r) + 5» a nf a S(r)] , $f 3) = -y/2 5^ J(r) , (2.5) 



A + G + — 

er 



where Latin indices run from 1 to 3 and /x runs from 1 to 4. 

Actually, there is some redundancy in this ansatz. A gauge transformation of the form 

A 5 = e *W oB (2.6) 

preserves the ansatz, but with new coefficient functions, which we indicate with a tilde, 
given by 

H + K = H + K, 

H — K — (P-ir)cos^ + 2 v / 2Jsin^, 
2V2J = -{H - Kjsinip + 2V2Jcosij , 
A-G = A-G, 

(a + G + — ) cos ip + 2V2F sin ip , 

2V2F = - ( A + G + simp + 2V2F cos ip , 
V er J 

S = S + F(1- cos,,) + ^ (A + G + I) sin xj) - ^ . (27) 

From the last of these equations, we see that S(r) can always be gauged away with a suitable 
choice of ip(r). We will henceforth assume that this has been done, so that S(r) vanishes 
identically. 

Requiring that the fields be nonsingular at the origin gives the boundary conditions 

A(0) = G(0) = H(0) = K(0) = . (2.8) 

The functions F and J can be nonzero at the origin. However, examination of the field 
equations, which we will display below, shows that nonsingular solutions must have 

F'(0) = J'(0) = , (2.9) 



2 In Rcf. [8], the third-sector fields were written in a somewhat different manner than here. The normaliza- 
tions in the present ansatz have been chosen so that the coefficient functions arc, nevertheless, the same 
as those that appear in that paper. 
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where a prime denotes differentiation with respect to r. 

To obtain the symmetry breaking that we want, the asymptotic value of the Higgs field 
must lie in the subgroup generated by the h a , giving the boundary conditions 

H(oo)=v, K(oo) = J(oo) = . (2.10) 

With this choice, the third-sector gauge fields are massive, so F(r) falls exponentially fast 
at large distance. The behavior of the other, massless, gauge fields depends on the magnetic 
charge. If the latter is a purely Abelian unit charge, then at large distance 

A(r) 1/er , G(r) < const/r 2 , F(r) ~ e~ evr/2 . (2.11) 

In Ref. [8] it was shown that in the BPS limit of vanishing scalar potential there is a 
solution given by 

A r = —J, ' 

sinn evr er 

H(r) = v cothetr , 

er 

v 1 



G(r) = K(r) = ( — ) L(r;b) 

V sinn evr er J 



F(r) = -J(r) = I v^M) > 

Vo cosh(et> r 1 1) 

S(r) = 0, (2.12) 
where b is any positive real number and 

L(r; b) = b — — . (2.13) 

v ; b + r coth(et;r/2) v ; 

This solution can be interpreted as being composed of two distinct fundamental 
monopoles. One is a massive monopole, with core radius ~ l/ev, whose magnetic charge 
has both non-Abelian and Abelian components. The other is a massless monopole that is 
manifested at the semiclassical level as a cloud of radius b whose magnetic charge cancels 
the non-Abelian part of the massive monopole's charge. This can be seen by computing the 
large distance behavior of the magnetic field. For 1/ev r b, both A{r) and G(r) fall 
as 1/r, and so 

^(D = ^ + °( 1 /r 3 ), (2-14) 

^( 2 ) = ^ + 0(l/r 3 ). (2.15) 

(The third-sector fields fall exponentially outside the massive monopole core and play no 
role here.) Outside the massless cloud, r>6, G(r) ~ —b/er 2 . As a result, Bf^ falls faster 
than 1/r 2 , while B?,^ is unchanged. Hence, the long-range magnetic field, and thus the total 
magnetic charge, have only first-sector components and are purely Abelian. 

One might think that the knowledge of the explicit analytic form of this monopole solution 
would make it an ideal case for constructing chromodyons. As we will see in the next section, 
this turns out not to be so. The difficulty arises from the fact that the energy of the BPS 
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monopole is independent of the cloud size b, so that a small perturbation can cause the 
cloud to expand without bound. To avoid this problem, we will add a potential term that 
effectively fixes the cloud size. With this term included, the BPS limit no longer applies, 
so we will have to solve the full set of second-order field equations. Because this cannot be 
done analytically, we will resort to numerical solution. 
Thus, let us add a potential of the form 

u 2 

V($) = -^-Tr$ 2 + a(Tr$ 2 ) 2 + Wr$ 4 (2.16) 

where /i, a, and b are constants. In order to obtain the desired symmetry breaking, b must be 
positive, while the requirement that the potential be bounded from below gives the condition 
4a + b > 0. At the minimum of the potential, 

Tr $ 2 = v 2 = . (2.17) 

4a + b v ' 

There is an SU(2) singlet Higgs scalar with mass 

m s = V2fi (2.18) 

and an SU(2) triplet with mass 

m t = y/2(l-c)ti, (2.19) 

where 

4a 

(2.20) 



4a + b 

The positivity of b implies that — oo < c < 1, with the two limits corresponding to m s <C m t 
and m t <C m s , respectively. 

Substitution of our spherically symmetric ansatz, Eq. (2.5), into the Euler-Lagrange 
equations, 

DjF ji = ie[$, £>*$], 
dV 

DjD 3 ® = — , (2.21) 

yields seven ordinary differential equations (ODEs), corresponding to the seven coefficient 
functions in our ansatz. Note that even though we can use the gauge freedom described 
by Eqs. (2.6) and (2.7) to make S(r) identically zero, there is still a corresponding ODE. 
However, while the other six ODEs are second order, this last is a constraint equation 
relating the coefficient functions and their first derivatives. If we set S = 0, this equation, 
which we will refer to as the S constraint, takes the form 







— e 



A + G + ^)f'-{a' + (?-^)f+\(H-K)J'-\(H'-K')J 

(2.22) 



1 f A + G + —\ 2 + -(H - K) 2 + 2J 2 + AF 2 



er J I 
while the other six ODEs become 
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A" = --A' + %A + 6eFF' + —A 2 + -H 2 

ry ry A ry ry 



+e 2 (A 3 + AH 2 + A J 2 + 5AF 2 - GF 2 - GJ 2 + 3HFJ + KFj) , 
G" = --G' + %G - QeFF' + —G 2 + -K 2 

ry ry£ ry ry 

+e 2 (G 3 + GK 2 + GJ 2 + 5GF 2 - AF 2 - AJ 2 - 3KFJ - HFJ) , 

F" = --F' + -(A + G)F' + e(A' + G')F + —(A + G)F+ —{H - K)J 
r 2 r 2r 

+e 2 ( AF 2 + A 2 + G 2 + -H 2 + -K 2 + 2J 2 + -HK - AG] F 
V 4 4 2 

e 2 

+ j(3AH - 3GK -GH + AK)J , 
H" = --H'+ %H - AeFJ' - 2eF'J + — AH 

ry ry £ ry 



+e 2 (3F 2 H + F 2 K + 2A 2 H + 6AFJ - 2FG.j) 
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V 



H(H 2 + 3K 2 - v 2 ) + AJ 2 (H - K) + c(4J 2 iT - 2iLFr 2 ) 



X" = —-if' + -^riT + AeFJ' + 2eF'J + — 

ry ry ry 

+e 2 (2G 2 if + 3F 2 if + F 2 # + 2AF J - 6GF j) 
+0l \k(K 2 + 3H 2 - v 2 ) + 4J 2 (if -H)- c(2H 2 K - 4J 2 ii) 

(■ 2.6 > 6 > • i 2c 2c 

J = — J' + -HF' - -KF' + eH'F - eK'F + —HF ifF 

r 2 2 r r 

2 2 

+— (A 2 + G 2 + 12F 2 - 2AG) J + — (3Aff - 3Gif + Aif - GH) F 



+ ,2 



2 

J 2 + (H - K) 2 -v 2 + c(2HK - AJ 2 )] J . (2.23) 



2_v > 2 

i> : 

These equations are not all independent. For example, the F" equation can be derived 
from the S constraint and the other five ODEs. The converse is not quite true, because 
there are solutions of the six second-order equations that do not satisfy the S constraint. 
However, if the S constraint holds at one value of r, the remaining ODEs imply that it holds 
for all r. In particular, the S constraint is satisfied for all r if the fields at spatial infinity 
obey Eqs. (2.10) and (2.11). 

In the BPS case, the exponential approach of the coefficient functions to their asymptotic 
behavior is governed by a single mass scale, ev. With the potential added, three different 
mass scales — ev, m s , and m t — come into play. To simplify our numerical simulations of the 
time evolution and to avoid the well-known stiffness problem, we want these characteristic 
lengths to be close to each other. To that end we set // = ev and c = 0.5, so that ev = m t = 
m s /V2. 

We choose to numerically solve the six ODEs in Eq. (2.23) and then check the solution 
against the S constraint. We use a MATLAB package - SBVP 1.0 [14] - to solve this 
boundary value problem. By using the collocation method, the SBVP numerical package 
can handle the singular terms in Eq. (2.23) with high accuracy near the origin, where 
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FIG. 1: Monopole solution in the non-BPS SO(5) gauge theory. The radial distance r is given in 
units of 1/ev, and the coefficient functions in units of v. 

r — > 0. Our numerical solution is shown in Fig. I. The S constraint has been checked to be 
automatically satisfied within the numerical error. 

It should be noted that this solution is not unique. Setting G = F= J = K = 
reduces the field equations to two coupled equations for A and H that are identical to 
those of the SU(2) theory. This then yields a solution that is simply an embedding of the 
SU(2) unit monopole via the subgroup generated by the h a . For our choice of parameters, 
the mass of this pure SU(2) solution is 1.287Mbps, where Mbps = 4:irv/e is the mass of 
the unit BPS monopole. By contrast, the non-embedding solution shown in Fig. 1 has a 
mass of 1.253M B ps- Note that the mass difference is numerically significant, well above the 
numerical errors. 

III. CONSTRUCTING A CHROMODYONIC CONFIGURATION 

In a U(l) dyonic soliton the electric charge results from a time-dependent phase of a com- 
plex field. In a similar fashion, a magnetically charged configuration with a time-dependent 
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orientation with respect to a non-Abelian group carries non-Abelian electric charge and is a 
chromodyon. In this section we will show how such chromodyons can be constructed from 
time-independent solutions such as those obtained in the previous section. We start with a 
static solution {(A)static(r )) ( ( ^ > ) static 

(r)}. As our first step, we excite one of the zero modes 
of the static solution and uniformly rotate its SU(2) orientation to obtain 

(A)l = R(t) (A)static# _1 (0 , 

($)i = R(t)($) static R-\t), (3.1) 

where 

R(t)= e ik3U0t (3.2) 

and the generator k% is defined in Eq. (2.3). It is critical to realize that this is not a gauge 
transformation, because the latter would have required that we also add a spatially constant 
A = —(u /e)k 3 . However, we do need a nonzero A in order to satisfy the Gauss's law 
constraint 

DjF j0 = ie[$, -0°$] . (3.3) 
Solving this equation, given (A)i an d ($)i and the boundary condition 

A)(oo) = 0, (3.4) 

yields a solution that we denote by (A ) n - This gives us a time-dependent configuration 
{(Ao) n , (A)i, (*)i}. 

It is often more convenient to work instead with the stationary configuration ob- 
tained by applying a gauge transformation with gauge function A = R" 1 . This gives us 
{(^o)in,(A) static j ($) static} where 

(A))m = — h + R~\A ) u R. (3.5) 

e 

It is not hard to see that we could have obtained this final configuration by directly solving 
Gauss's law with (A) static and ($) static given and Aq obeying a different boundary condition, 

A (oo) = ^%. (3.6) 
e 

We denote the solution of this equation as (A))static- 

Even when working in this static gauge, it is convenient to describe the extra energy 
associated with this configuration in rotational terms. The relevant "moment of inertia" is 
given by the spatial integral of the sum of the squares of the field components that undergo 
the phase rotation. Since the characteristic length scale is ~ {ev while the natural scale 
of the fields is v, we have 

^ch~^ 2 ~f-)% 2 c 2 ~-f^) 2 . (3.7) 

\evj e\evj 

Similarly, the SU(2) electric charge has a magnitude 

q E ~ elu ~ - (3.8) 



e \ev 
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with the extra factor of e arising because the electric charge is e times the momentum 
conjugate to the phase rotation. 

The configuration {(A ) st3itic , (A)static, ($)static} has both a U(l) magnetic charge and an 
SU(2) electric charge, and so is a chromodyon. The question that we need to address is 
whether it is a solution of the field equations. It is easy to see that it cannot be, because 
the rotation in group space induces a deformation of the field profiles, just as the spatial 
rotation of a solid object induces a deformation of its shape. However, this deformation is 3 
of order Uq, and so should require only a small modification of the configuration if ojq is 
sufficiently small. If this is the only correction needed, then the theory does indeed have 
a stable static chromodyon solution. On the other hand, it may be that there is no such 
static solution. We will test for this possibility by taking {(A ) s tatic, (A)static, ($)static} as an 
initial condition and then letting the fields evolve in time according to the field equations. 
If the deformation resulting from the rotation in SU(2) space is the only impediment to its 
being a static solution, the fields should oscillate, with an initial amplitude proportional to 
lVq, and eventually settle down in the true static solution. 

In order to do this, we need to determine (A))static- Thus, our immediate task is to solve 
Eq. (3.3) subject to the boundary condition that A (oo) = (uj /e)k 3 . The fields (A)static 
and ($)static are both spherically symmetric. However, the boundary condition on A breaks 
this symmetry, so we can only assume that A has an axial symmetry. The most general 
ansatz for A is then 



A)(l) 


= r a u{p, z) , 


A l(i) = w{p,z) , 


^0(2) 


= r a b(p,z), 


A)(2) = Q{Pi z ) 5 


\a 

A 0(3) 


= V2e ab f b q(p,z) , 


^0(3) = > 


A 4 


= -V2t(p, z), 





(3.10) 

where a and b are either 1 or 2 and p = \Jx 2 + y 2 . The boundary conditions at spatial 
infinity require that Q approach u /e and that the other coefficient functions tend to zero. 

Substituting this ansatz into Gauss's law, Eq. (3.3), yields the set of second-order partial 
differential equations that are displayed in Appendix A. These rather complicated equations 
experience a remarkable simplification if (A) static and ($) sta tic are taken to be the BPS 
solution of Eq. (2.12). In this case, they can be satisfied by setting all the coefficient 
functions except Q to zero, taking Q(p, z) = Q(r), and requiring that 

^ + eG(r)Q(r) = 0, (3.11) 
where G(r) is the second-sector gauge field function given in Eq. (2.12). The solution to 



3 This can be seen easily from the field equations 

D Q F° l + DjF' ji = ie[$, D l §] , 

D a D <I> - DjDj<f> = 0, (3.9) 

where the non-static terms, D F 0t and D D <I>, are of order 0(oJq). If these terms are omitted, the 
equations reduce to the static equations satisfied by (Aj) stat i c and (^static- 
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Eq. (3.11) is 4 

n( ) UJ ° b = fu \ rcoth(evr/2) 

W[T) eL(r,b) db V e / b + r coth(et;r/2) ' l ' j 

Unfortunately, this simple solution turns out not to be useful. To see this, recall that 
for small velocities the dominant time dependence of a soliton arises entirely through exci- 
tation of its zero modes, whose dynamics is governed by the moduli space Lagrangian. For 
the generic case, with a nonzero scalar field potential, there are seven zero modes about 
{(^4j)static, (^static}- Four of these — three translation modes and one U(l) phase mode 
- are irrelevant for our purposes. The remaining three are SU(2) orientation modes, one 
of which has been excited by the transformation in Eq. (3.1). Within the moduli space 
approximation (MSA) there would be uniform motion in the corresponding collective coor- 
dinate. In the gauge where A (oo) = 0, the soliton would rotate uniformly in SU(2) space, 
as described by Eq. (3.1). 

In the BPS case there is an additional zero mode, corresponding to the freedom to vary 
the cloud radius b. The moduli space Lagrangian governing the eight zero modes is given 
by 



L M s = ImX 2 + 



2 2M 




d 2 + sin 2 a f3 2 + (j + cos a $) 2 \ , (3.13) 



where M is the BPS monopole mass, X is the location of the center of the system, x is the 
U(l) phase, and a, (3, and 7 are SU(2) Euler angles. We are interested in the part within the 
curly braces, L^ s , that describes the zero modes corresponding to the non-Abelian cloud 
size and the SU(2) orientation of the non-Abelian cloud. By transforming to coordinates 

n rr a 3 — 7 

x 1 = 2vosm— cos — - — , 

Z Z 

rr a 3 + 7 
x 2 = 2vocos— cos — - — , 

n rr • a ■ P ~ 7 
x 3 = 2vosm— sin — - — , 

1— Ol 3 + T 

x A = 2Vfecos-sin- -, (3.14) 

we see that this is actually the Lagrangian for a free particle in four-dimensional Euclidean 

space, 

-^ms = ~*~ 2^ 2^ • (3.15) 

The solutions of this Lagrangian are uniform straight-line motion. Without loss of gen- 
erality, we focus on solutions in the X\-x 2 plane. The SU(2) rotating configurations we are 
studying then correspond to taking initial values a = 3 = 7 = 0, b = bo, a = u>o, and 
ft = 7 = b = 0. This leads to 

Xl(t) = v t, 

X2{t) = Po, (3.16) 



This result was also obtained, by a different method, in Ref. [15]. 
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or, equivalently, 



m = bo+b (^) , 

a(t) = 2 tan" 1 ^. (3.17) 

We see that the SU(2) phase does not even go through a full rotation, so this hardly a good 
approximation to a uniformly rotating configuration of fixed color-electric charge. This is 
clearly attributable to the fact that the cloud size can grow without limit. 5 

To avoid this difficulty, we turn to the case with a nonzero potential, where the b mode 
is no longer a zero mode and the MSA predicts uniform rotation of the SU(2) orientation. 
Because analytic results are no longer possible, we must resort to numerical solution of the 
field equations, both to obtain the static monopole solution, as described in the previous 
section, and to find the A that solves the Gauss's law constraint, the topic to which we now 
turn. 

The ansatz for Aq was given in Eq. (3.10), and the coupled field equations that follow 
from this ansatz are given in Appendix A. The outer boundary conditions are found by 
noticing that well outside the core [i.e, when r is much greater than (ev m' 1 , and m^ 1 } 
the third-sector components are all exponentially small and we have 

u(p,z) 

w(p,z) 

b{p, z) 

Q(p,z) 

Here c\ and c 2 are free constants, to be determined from the numerical simulation, while d\ 
and d,2 can be derived in terms of these by analysis of the asymptotic expansion. We use the 
successive over-relaxation (SOR) method with red-black ordering as our numerical method. 
Our results are shown in Fig. 2. 

In contrast with the BPS case, we see that all of the coefficient functions are nonzero, 
although Q remains dominant, and that all of these functions, including Q, have only axial 
symmetry, with separate dependence on p and z. 

IV. EVOLVING THE CHROMODYON 

In the previous two sections we obtained a static monopole solution and then determined 
the A that is required by Gauss's law when this solution rotates uniformly in SU(2) space. 
We now take this configuration as the initial condition and let the system evolve as dictated 
by the equations of motion. We work in the gauge where the uniform rotation has been 
gauged away, so that the initial configuration {(A ) sta tic, (A)static, ($)static} would be a static 



5 Equation (3.17) implies that b increases linearly with time. In actual fact, the MSA breaks down when 
when b approaches the speed of light [16]. 



dpz 



C\Z 

dipz 



C2 d 2 d\z 2 



oJo H h -J + , 



(3.18) 
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FIG. 2: Numerical solution for (^4o)static in the non-BPS SO(5) theory, with ojq = 0.1 ev. The 
coordinates z and p are given in units of 1/ev and the component fields u, w, b, Q, q, and t, defined 
in Eq. (3.10), in units of v. 

solution if the MSA were exact. In this gauge, any time dependence arises from corrections 
to the MSA. 

To proceed, we need to specify u . It cannot be too big (e.g., so large that the energy 
arising from the gauge rotation is comparable to the monopole mass) if the original configu- 
ration is to be even an approximate solution. A more stringent condition is suggested by the 
existence of the embedded pure SU(2) monopole described at the end of Sec. II. In order 
to make sure that our configuration does not evolve toward this other monopole solution, 
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we want the energy associated with the gauge rotation to be less than the mass difference 
between the two types of static monopoles. On the other hand, taking ujq to be too small 
will impose increased computational burdens, because we will have to simulate the evolution 
for a much longer time in order to see any effect. 

We choose u = 0.04 ev. The gauge rotation energy corresponding to this is of order 
(uj /ev ) 2 Mbps ~ 10~ 3 Mbps, smaller than the mass difference between the two monopole 
solutions. Because the A obtained in the previous section is linearly proportional to u , 
the initial data can be obtained by a simple rescaling of the solution shown in Fig. 2. 

The Euler-Lagrange equations consist of the evolution equations 

D F 0i + DjF ji = ie[$, D%] , (4.1) 
D D°$ + DjEPQ = , (4.2) 

and the Gauss's law constraint 



DjF j0 = ie[<E», _D°$] . (4.3) 

We will use the so-called free evolution scheme to simulate this constrained system. In this 
scheme, we numerically solve the evolution equations, and use the constraints to monitor 
the accuracy of the evolution. It is well known that direct numerical implementation can 
have problems with numerical instability. To avoid this, we choose to use the technique in 
Ref. [17] to first rewrite the evolution equations in hyperbolic form. To do this, we take 
a covariant time derivative of Eq. (4.1) and subtract a covariant spatial derivative of the 
Gauss's law constraint, obtaining 

D°D F 0i + D°DjF ji - D'DjF^ = ieD°[$, - ie D°$] . (4.4) 

Switching the two covariant derivatives on the left-hand side gives 

D°D F 0i + DjD F ij + DjDiFjo + 2ie[F ij , F j0 \ = ie D°[<Z>, D% - ie £>*[$, _D°$] . (4.5) 

Next, we take a covariant spatial derivative of the Bianchi identity, 

D Q Fij + D t F j0 + DjFoi = , (4.6) 

to give 

DjD Q Fij + DjDiFjo = DjDjF i0 . (4.7) 
We substitute this into Eq. (4.5) and obtain 

D°D F i0 + D j DjF i0 + 2ie[F ij ,F j0 ] = ie - ie £>*[$, D°$] . (4.8) 

The full set of equations now consists of the definition of F^ u and two wave equations, 
Eqs. (4.2) and (4.8). Equation (4.1) becomes a second constraint. 

One advantage of this hyperbolic formulation is that it does not depend on the gauge 
condition. Therefore, by assuming an appropriate time- dependent gauge transformation, we 
can choose A to be anything we want. In particular, we will set 

A (r,t>0) = (A (r)) static , (4.9) 
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FIG. 3: (a) Using three slabs in the three-dimensional Cartesian coordinates, (b) Axial symmetry 
is used to update grid points on the off-center slabs. 

where (A ) s tatic is the solution shown in Fig. 2. The chromodyon configuration constructed 
previously is an approximate stationary solution with this choice of A . If it were exact, 
Eq. (4.9) would correspond to choosing a gauge in which the SU(2) rotation was gauged away. 
Any time dependence that we observe will correspond to corrections to this approximation. 

Once A has been fixed in this manner, the time-dependent variables in this formulation 
are A i: $, and Ei = F i0 . Their initial values at t — are chosen to be (A)static, (^static, 
and Dj(Ao) static- The time derivatives of Ai and $ at t = are set equal to zero, while 
doEi(t = 0) is obtained from the constraint Eq. (4.1). Note that Ei(t = 0) is 0{ujq) and 
d Ei(t = 0) is 0(ujI). 

The system has axial symmetry at t — and this symmetry will be preserved by the evolu- 
tion. Usually, a two-dimension grid structure in the p-z plane is used to discretize an axially 
symmetric system. However, using cylindrical coordinates to simulate time-dependent sys- 
tems can easily cause numerical instabilities [18]. In our implementation, we use the method 
proposed in Ref. [19]. The idea is to discretize the system using three slabs in the three- 
dimension Cartesian coordinates, as shown in Fig. 3a. Each slab is a two-dimensional grid 
structure with step size h. The central slab lies in the x-z plane. The other two slabs are 
obtained by shifting the central one by +h and —h along the y axis, respectively. At each 
grid point, any so(5)-valued element is represented by ten real numbers. For each time step, 
the two wave equations are solved first on the central slab using the finite difference method. 
Then the axial symmetry is used to update the grid points on the other two slabs. To be 
precise, we plot the grid structure in Fig. 3b. The axial symmetry tells us that the system 
will be invariant if we rotate it through the same angle both in real spatial space and in 
isospin space. For the $ field, this means that 

$(r D , t) = e- ihs9 e- ikie <5>{r E , t)e iki9 e ihid , (4.10) 

where 9 is the angle between r E and r D . For the A; t field, it means that 

Mr D ,t) = e-™e-™R ij (9)A j (T E ,t)e ikl > e ih * e , (4.11) 
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with Rij(0) being the matrix corresponding to a spatial rotation by angle 9 about the z-axis. 
Interpolation is used to calculate $(r^,t) and Aj(i\E,i) from the values at the neighboring 
grid points, A, B, and C. 

In our simulation, the grid with three slabs covers a space of size 60 (in units of 1/ev). 
The open boundary condition is used because of its simplicity. The radius of the monopole 
is roughly 1 and we focus on studying a region with a radius of 10. This implies that we 
can only simulate our system up to t = 50. 

Because the system is approximately stationary at t — 0, the time-dependent parts are 
very small. In the hyperbolic formulation, we have four fields, A , A iy Ei = F i0 , and <3>. 
Since we have already prescribed Aq to be time-independent by Eq. (4.9), only Aj, Ei, and 
$ have time- dependent parts. In our implementation, we separate out the time-dependent 
parts via 

Ai(r,t) = Ai(r)+ A(r,t), 

$(r,t) = $(r) + $(r,t), (4.12) 

where Aj(r, t = 0) = E^r, t = 0) = $(r, t = 0) = 0. During the evolution, we check how well 
these time-dependent parts satisfy the Gauss's law constraint. Theoretically, this constraint 
should be satisfied at t — because of the way we construct the initial gauge-rotating 
system, and should remain satisfied for all t. Numerically, however, it is only satisfied up to 
some finite accuracy. To check how well the Gauss's law constraint is satisfied, we write the 
right-hand side of the constraint as a sum of nine terms, and then treat each of these terms 
as a vector in a ten-dimensional space. 6 We then define the sum of all these vectors (which 
theoretically should vanish) to be the defect. This defect should remain small as long as our 
numerical scheme is stable. If the scheme is not stable, the defect will grow exponentially. 
We calculate the ratio between the maximum norm of the defect and the maximum norm 
among all of its component vectors. (For both norms we take the maximum over the entire 
computational domain.) We use this ratio to measure the accuracy to which the Gauss's 
law constraint is satisfied. At t — this ratio is 0.5 x 10~ 2 , while during the evolution the 
ratio for the time-dependent parts never exceeds 3 x 10~ 2 . More details about this can be 
found in Appendix B. 

V. GLOBAL GAUGE ROTATION SLOW-DOWN 

Our initial configuration satisfies the equations of motion to first order in iv . If this 
first-order approximation were exact, the solution would be time-independent, so the time- 
dependent parts in Eq. (4.12) tell us how the real evolution deviates from this first-order 
approximation. By analyzing the result of our simulation, we find that we can best fit 
this deviation in terms of a global gauge rotation. Recall that our initial configuration 
corresponds to a monopole rotating in SU(2) space about the k 3 axis with an angular velocity 
ujq. We made it static by going to a gauge with A (oo) = (c<Jo/e)A;3, effectively transforming 
to a rotating frame. By fixing A (r,t) as in Eq. (4.9) we stay within that rotating frame 



This is because every term in the Gauss's law constraint is in the ten-dimensional SO (5) adjoint repre- 
sentation. 
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during the evolution. Any slowing of the gauge rotation would appear in this frame as a 
rotation about the fc 3 axis, but in the opposite direction. Thus, it would correspond to a 
global gauge rotation generated by 

A s (0) = e -^« (5.1) 
with positive 6{t). The effect would be as if the initial angular velocity c^o were replaced by 

d9 



dt ' 



(5.2) 



leading to a configuration with smaller color charge and a smaller energy. 

The fields Ai and $ can be expanded in SO (5) components, as in Eq. (2.4). The first- 
sector components are unchanged by the rotation generated by A s (0). The second-sector 
components transform as 





( 


p2 






V 



/ cos 6* 



sintf 0\ ( Pl 2) \ 
sin 9 cos 9 
1 j 



(2) 
p2 
M2) 

V J 



(5.3) 



while the third-sector components decompose into two doublets transforming according to 
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Let us define 



A$(r,t;0) 
AA,(r,^) 



A^^A; 1 ^)-^!-) -$(r,f), 



A^A^A; 1 ^)-^) 



(5.5) 



where $(r,t) and Aj(r,t) are from our numerical simulation. If the time evolution of our 
configuration were completely due to a slowing of the global gauge rotation, then for any 
given time t there would be a single 9(t) that would make A$(r, t; 9) and AAj(r, t; 9) both 
vanish for all values of r. To see how close we are to this situation, we can extract a value 
for 9 by several different methods and then compare these values. First, we obtain 9 from 
the second-sector components of $ by minimizing 



Tr 



A$ (2) (r,f;0) 



(5.6) 



at various points. In Fig. 4a we show the 9 obtained in this manner for a series of points 
along the x-axis. As can be seen, the #'s thus obtained are only weakly position-dependent. 
We can also define a spatially averaged 9 by finding the value that minimizes quantities such 

as 

Af(t, 9)= J fJ^ R {Tr [AA(r, t; 9)f + Tr [A$(r, t; 9)} 2 } dxdydz , (5.7) 
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(c) (d) 

FIG. 4: (a) Fitted values of 6 at grid points along the positive x-axis. (b) Values of d9/dt obtained 
by fitting with various sets of component fields. Circles indicate fits using $, squares fits using 
Ai, the solid line fits using A^ 3 \, and the dashed line fits using A^y (c) Comparison of results 
using different grid sizes. The long dashed curve, the dashed curve, and the solid line curve are 
calculated using Ah = 0.18(ev) _1 , Ah = 0.12(ev) _1 , and Ah = 0.08(eu) — 1 grid sizes, respectively, 
(d) Plot of rj(t), defined by Eq. (5.8), which indicates the fraction of the time dependence that 
cannot be accounted for by the global gauge rotation. In all four figures distances and times are 
given in units of 1/ev and ujq = O.OAev. 

where R e is the size of the physical region of interest. 7 In Fig. 4b, we compare the results for 
dO/dt that are obtained by restricting M in several different ways: using only $ or only A^, 
or using just the second-sector or just the third-sector components of A { . [The first-sector 



7 The numbers we present are obtained using R e = 6(et>) 1 , but these results are not very sensitive to the 
exact value of R e . 
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FIG. 5: Decomposition of d9/dt (solid line) into linear and oscillating components. Time is in 
units of 1/ev, and ujq = 0.04et> . 



components are invariant under A s (6), and so cannot affect the fitting of 8.} We see that all 
of these methods give essentially the same results, again consistent with the interpretation 
in terms of a spatially uniform global gauge rotation. In order to indicate the convergence 
of our simulations, in Fig. 4c we show three curves for uq = 0.04et> , using successively finer 
grid structures. 

We can also ask how much of the time-dependence can be accounted for by this uniform 
rotation. To this end, we divide the residual norm after fitting 9 by the norm of all the 
component fields in the time- dependent parts, and define 



V (t) = 



Af(t,0) 



(5.8) 



We plot i] for our simulation in Fig. 4d. We see that, although initially the global gauge 
rotation accounts for only a small part of the time dependence, at large times it is clearly 
the dominant component. 

We see from the data in Fig. 4 that d6/dt increases with time (although not uniformly), 
with a corresponding decrease in cu cff . This can be interpreted as the sum of two effects, as 
shown in Fig. 5. One is an overall oscillation pattern that appears to be a transient effect 
caused by the relaxation of the system after the initial excitation. The other effect is a 
linearly increasing d9/dt. By the end of our simulation, at t = 52(et>)~ 1 , about 5% of the 
initial angular velocity has been lost. We do not see any indication that this slowing down 
process will stop, and expect that u c g would eventually tend to zero if the simulation could 
be carried out for long enough. 

The energy that the chromodyon loses through the slowing of the gauge rotation must be 
carried away by radiation of the massless gauge fields in the unbroken non-Abelian subgroup. 
(Indeed, much of the nonrotational contribution to i] at early times can presumably be 
attributed to the creation of the radiation field.) It is of interest to know how this radiation 
depends on uj . We can use the slope of the constant part in Fig. 5 to approximate d 2 9/dt 2 for 
uj = 0.04et> . In Table I we show the results for this as well as the corresponding results for 
simulations with two other values of u . Recalling from Eq. (3.7) that the energy associated 
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d 2 6/dt 2 


(ev/u)l)d 2 6/dt 2 


0.01 ev 


8.0 x 10~ 7 (ev) 2 


0.80 


0.02 ev 


6.2 x 10~ 6 (ev) 2 


0.78 


0.04 ev 


4.6 x 10~ 5 (ev) 2 


0.72 



TABLE I: Dependence of the deceleration of the initial global gauge rotation on the initial angular 
velocity loq. 



with the phase rotation is proportional to or, we see that this data is consistent with 

dE ch du cS 4 

-w (5 - 9) 



We can understand this dependence by considering the energy flux carried by the radiation 
of the massless non-Abelian gauge fields in the unbroken subgroup. This should be given by 
the analogue of the electromagnetic Poynting vector, 

Ti ~ e ijk Tr EjB k , (5.10) 

where the hats indicate the 0(l/r) radiation components of the field strengths. Because 
the initial chromodyon configurations that we constructed satisfied the static field equations 
to first order in cu , these radiation fields must be each at least second order in uo, so that 
Ti ~ ufi. 



VI. CONCLUDING REMARKS 

Magnetic monopoles can be promoted to dyons by time-dependent excitation of their U(l) 
global gauge zero modes. In this paper we have addressed the question of whether monopoles 
in theories with non-Abelian unbroken symmetries can be promoted to chromodyons - 
monopoles with non-Abelian electric charge — by a similar excitation of their non-Abelian 
global gauge zero modes. It has long been known that the answer is negative if the magnetic 
charge has a non-Abelian component, because there are then topological obstructions that 
preclude the existence of a chromodyon. However, there are also monopoles with purely 
Abelian asymptotic magnetic charge, for which there is no such obstruction, that could 
potentially have chromodyonic counterparts. We have examined one such case here, using a 
constructive approach. We started with a configuration with a globally rotating non-Abelian 
phase, and thus a nonzero chromo-electric charge, and then numerically evolved it to see 
whether it would settle down in a stable static solution. In our simulations we found instead 
that the effective rate of gauge rotation slows down, so that the chromodyon continually loses 
energy and chromo-electric charge. Although we were not able to continue the simulation 
until this charge was completely lost, every indication suggests that this would be the final 
state of the system. 

It is instructive to compare our results with those that would have been obtained by 
applying our methods in the theory with SU(2) broken to U(l), where we know that there 
is a dyon with Abelian electric charge. Because there is no analogue of the cloud radius 
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zero mode, with its associated complications, we can work in the BPS limit, where analytic 
expressions are available. 

Our approach would start with the static monopole solution 



with 



$ a = f a H(r), 



Mr) = — r , 

sinh evr er^ 

H(r) = fcothetr . 

er 



(6.1) 



(6.2) 



Applying a global U(l) phase rotation and then gauge transforming back to a static gauge 
would yield 

(A$) static = r a Q(r) , 



where 



Q(r) 



U! 



ev 



v coth evr — 



er 



From the 1/r term in this expression, we see that the asymptotic electric field is 



E a = f f — 



where 



The energy of this configuration is 



1 /^o\ 



e Vet;/ 



E 



4ttv 



1 + 



e 2 q 2 E 



(6.3) 
(6.4) 

(6.5) 
(6.6) 

(6.7) 



where the first term represents the mass of the original monopole and the second is the 
additional energy due to the phase rotation. 

As in the chromodyon case, this initial configuration is only an approximate solution of 
the equations of motion. The exact dyon solution with charge q E is given by [20, 21] 



A(r) 



1 



sinh ev'r er 
H{r) = cosh 7 



v coth ev'r 

er 



Q(r) = sinh 7 



v' coth ev'r 

er 



(6.8) 



where 7 is determined by the ratio of electric and magnetic charges and is given by 

sinh 7 = eqE (6.9) 
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and v' = v / cosh 7. It has an energy 



E = —y/l + (eq E )* (6.10) 
e v 

and corresponds to phase rotating with an angular velocity 

' • v, e2vqE faii\ 

uj = ev smh7 = = . (o.ll) 

Thus, our construction would start with a configuration that has a core radius that is 
a factor \Jl + (eg^) 2 smaller than that of the exact solution, and an angular velocity that 
is larger by the same factor. (The smaller core radius produces an decrease in the phase 
rotation moment of inertia that exactly compensates for the increase in angular velocity, 
thus yielding the same electric charge.) The energy of this initial configuration exceeds that 
of the exact dyon solution by an amount of order q E . 

It is easy to see what will happen if the initial configuration of Eqs. (6.2) and (6.4) is 
allowed to evolve. Because radiation of the massive charged gauge field is energetically 
suppressed, the electric charge of the dyonic configuration will be conserved. Hence, as 
the initial system relaxes it will tend toward the static dyon solution of the same charge. 
It will shed energy by radiating massless photons, but the amount of this energy loss is 
constrained by the fact that exact dyon mass places a lower bound on the energy. As the 
dyon radiates its phase rotation will slow down and its core will expand. For small electric 
charge (i.e., eq E <C 1), this slowing and expansion will both be small, and the system will 
quickly approach its final state. In particular, the slowing of the phase rotation will be far 
less than that which we found in our non-Abelian simulation. 

Thus, our numerical simulations provide strong evidence against the existence of static 
chromodyons in a theory with SO(5) broken to SU(2)xU(l). Because it is hard to see 
how enlarging the unbroken symmetry to a different non-Abelian group would stabilize the 
chromodyon, we expect that similar results would hold for other choices of gauge group 
and symmetry breaking. Of course, numerical simulations cannot provide a rigorous proof. 
Even apart from issues related to numerical accuracy, there is always the possibility that 
the specific choice of initial configuration played a crucial role. For example, it is logically 
conceivable (although we think it quite implausible) that there is some special choice or 
range of u that would have led to a stable chromodyon. Another possibility is that there 
are chromodyon solutions, but that these exist only for some minimum value of the chromo- 
electric charge. In this case, the solutions would not be continuously related to the purely 
magnetic monopole, and so might not be found by our method. Although we cannot exclude 
this possibility, it seems to us to be rather unlikely. Hence, subject to these caveats, we 
conclude that static chromodyon solutions do not exist. 
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APPENDIX A: GAUSS'S LAW EQUATION WITH AXIAL SYMMETRY 

As discussed in Sec. Ill, when we gauge rotate the spherically symmetric static non-BPS 
monopole solution, the resulting Gauss's law equation has only axial symmetry and consists 
of six coupled partial differential equations for the six coefficient functions u, w, b, Q, q, 
and t appearing in the A ansatz of Eq. (3.10). We write out the detailed forms of these 
equations below. Here, A, H, G, K, F, and J are the functions, defined by Eq. (2.5), that 
specify the static non-BPS monopole solution and that are shown in Fig. 1. To simplify the 
equations, we have set e = 1 throughout; the explicit factors of e can be recovered by simple 
dimensional analysis. 

The two equations corresponding to first-sector components are 



j ii o z 

d pp u + -d p u - — + d zz u - 2^Ad z w + 2-Ad p w - AFd z q + AFd p t 

+2-(d r F)t - 2-(d r F)q - Q-AFt + 2-FGt - ^KJt + ^HJt 

ry ry ry ry ry ry 



+2-AFq - 2-FGq + -KJq + 3-HJq - 2 P -A 2 u - ^A 2 w + ^H' z w 
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(Al) 
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P 2 u2 



. z 



(A2) 



The two second-sector equations can be obtained from these simply by making the sub- 
stitutions 



u — > b , w — > Q , 
A^G, H^K. 



(A3) 



Finally, the third-sector equations are 

d PP q + -d p q - 4 + d zz q + -Ad z t - ^Gd z t + Fd z b + Fd z u - -Ad p t + -Gd p t 



-Fd p Q - Fd p w - 2 P -{d r F)Q - \ p -{d r F)w + \-{d r F)b + \-(d r F)u 
r 2 r 2 r 2 r 



r 



-AFQ - -FGQ + -HJQ + -JKQ - -AFw + -FGw - -HJw - -JKw 
2 V 2 V 4 V 4 V 2 2 4 4. 



--AFb + -FGb - \ujb - %JKb + -AFu - -FGu + ^rHJu + \ JKu 
22442244 



z 

+ - 

r 

- l -A 2 q - X -G 2 q - - A H 2 q - hlKq - l -K 2 q - AF 2 q - 2J 2 q - l -Aq - -Gq = , (A4) 
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APPENDIX B: MONITORING THE GAUSS'S LAW CONSTRAINT 



Analytically, the Gauss's law constraint should remain satisfied throughout the evolution 
if it is satisfied by the initial data. In our numerical calculation, we keep monitoring how 
well the constraint is satisfied and use this as a way to check the accuracy of our numerical 
methods. 

The Gauss's law constraint is 



DjF j0 



ie[$,D°$] . 



(Bl) 



We can use Eq. (4.12) to expand this into time- dependent and time- independent parts. 
Focusing on the time-dependent part, we have the requirement that 

= djEj + ie[Aj, Ej] + ie[Aj, Ej\ + ie[Aj, Ej] - ie[<l, d i>] - ie[4>, <9 $] 

-ie 2 [§,i[AQ,§}} -ie 2 [$,i[A ,$]] - ie 2 [$>,i[A ,§}} . (B2) 

To demonstrate that the Gauss's law constraint is satisfied in our numerical calculation, 
we calculate these nine terms using our numerical results and then define their sum to be 
the defect. Since the defect is an so(5) element, we can define its norm to be the square root 
of the trace of its square. In Figs. 6 and 7 we plot on the x-z plane the norm of the defect 
and, for comparison, the norm of the first term, both at t — 34.2 (ev)" 1 . As we can see, 
the defect is about two orders of magnitude smaller than the first term, indicating that the 
Gauss's law constraint is satisfied very well. We define the defect level to be the ratio of the 
maximum norm of the defect and the maximum norm among all terms. (Both maxima are 
found over the entire computation domain.) For example, at t — 34.2(et> ) _1 , the maximum 
norm of the defect is about 8 x 10~ 6 , while the maximum norm of the first term (which is 
larger than that of any of the other terms) is about 4 x 10~ 4 . This give a defect level of 
about 0.02. 



[1] A. Abouelsaood, Nucl. Phys. B 226, 309 (1983). 

[2] P. C. Nelson and A. Manohar, Phys. Rev. Lett. 50, 943 (1983). 

[3] A. P. Balachandran, G. Marmo, N. Mukunda, J. S. Nilsson, E. C. G. Sudarshan and F. Zac- 
caria, Phys. Rev. Lett. 50, 1553 (1983). 



24 



FIG. 6: The norm of the defect of the Gauss's law constraint at t = 34.2(et>) _1 . Since the system 
has an axial symmetry along the z axis, we plot the results on the x-z plane. Distances are given 
in units of 1/ev, and the defect is in units of v. 




FIG. 7: The norm of the first term in Eq. (B2) at t = 34.2(e?;)~ 1 . Since the system has an axial 
symmetry along the z axis, we plot the results on the x-z plane. Distances are given in units of 
1/ev, and the defect is in units of v. 

[4] A. P. Balachandran, G. Marmo, N. Mukunda, J. S. Nilsson, E. C. G. Sudarshan and F. Zac- 

caria, Phys. Rev. D 29, 2919 (1984). 
[5] A. P. Balachandran, G. Marmo, N. Mukunda, J. S. Nilsson, E. C. G. Sudarshan and F. Zac- 

caria, Phys. Rev. D 29, 2936 (1984). 
[6] P. A. Horvathy and J. H. Rawnsley, J. Math. Phys. 27, 982 (1986). 
[7] P. A. Horvathy and J. H. Rawnsley, Phys. Rev. D 32, 968 (1985). 
[8] E. J. Weinberg, Phys. Lett. B 119, 151 (1982). 



25 



[9] B. Julia and A. Zee, Phys. Rev. D 11, 2227 (1975). 

[10] G. 't Hooft, Nucl. Phys. B 79, 276 (1974). 

[11] A. M. Polyakov, JETP Lett. 20, 194 (1974) [Pisma Zh. Eksp. Teor. Fiz. 20, 430 (1974)]. 

[12] P. C. Nelson, Phys. Rev. Lett. 50, 939 (1983). 

[13] P. C. Nelson and S. Coleman, Nucl. Phys. B 237, 1 (1984). 

[14] W. Auzinger, G. Kneisl, O. Koch, and E. B. Weinmuller, Numerical Algorithms 33, 27 (2003). 

[15] K. Lee, E. J. Weinberg and P. Yi, Phys. Rev. D 54, 6351 (1996) 

[16] X. Chen, H. Guo and E. J. Weinberg, Phys. Rev. D 64, 125004 (2001) 

[17] A. Abrahams, A. Anderson, Y. Choquet-Bruhat and J. W. . York, Phys. Rev. Lett. 75, 3377 
(1995) 

[18] L. Pogosian and T. Vachaspati, Phys. Rev. D 62, 105005 (2000) 

[19] M. Alcubierre, S. Brandt, B. Brugmann, D. Holz, E. Seidel, R. Takahashi and J. Thornburg, 

Int. J. Mod. Phys. D 10, 273 (2001) 

[20] M. K. Prasad and C. M. Sommerfield, Phys. Rev. Lett. 35, 760 (1975). 

[21] E. B. Bogomolny, Sov. J. Nucl. Phys. 24, 449 (1976) [Yad. Fiz. 24, 861 (1976)]. 



26 



